IPMU 09-0052 



Phenomenological Aspects of Horava-Lifshitz Cosmology 

Shinji Mukohyama*-"^ Kazunori Nakayama*^*-' , 
Fuminobu Takahashi^"\ Shuichiro Yokoyama*^'^-* 
° Institute for the Physics and Mathematics of the Universe, 
University of Tokyo, Chiba 277-8568, Japan 
^ Institute for Cosmic Ray Research, 
University of Tokyo, Kashiwa 277-8582, Japan, 
Department of Physics and Astrophysics, 
Nagoya University, Aichi 4^4-8502, Japan 
(Dated: May 7, 2009) 

Abstract 

We show that, assuming the dispersion relation proposed recently by Hofava in the context of 
quantum gravity, radiation energy density exhibits a peculiar dependence on the scale factor; the 
radiation energy density decreases proportional to a^^. This simple scaling can have an impact on 
cosmology. As an example, we show that the resultant baryon asymmetry as well as the stochastic 
gravity waves can be enhanced. We also discuss current observational constraint on the dispersion 
relation. 
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I. INTRODUCTION 



Recently a new class of quantum gravity was proposed by Hofava motivated by the 
solid state physics. This theory is fundamentally non-relativistic and, in the ultraviolet 
(UV), exhibits the Lifshitz scale invariance 



t hH. 



X 



(1) 



with dynamical critical exponent z = 3. It is this anisotropic rescaling that makes Horava's 
theory power-counting renormalizable. The relativistic scaling with z = 1 is restored in the 
infrared (IR) due to deformation by relevant operators, and the Lorentz symmetry emerges 
as an accidental symmetry. 

The z = 3 scaling implies that, in the UV, a physical degree of freedom should have a 
dispersion relation of the form 

where M is a characteristic mass scale. Since the Planck mass is an emergent quantity in 
Horava's theory, M does not have to be the same order of magnitude as the Planck mass. 
Moreover, M can differ for different physical degrees of freedom: graviton, photon, quarks, 
leptons, scalar fields, etc. The fundamental symmetry of Horava's theory is invariance under 
foliation-preserving diffeomorphism: 
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x'{x^,t), t^i{t). 
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Any value of M is consistent with this symmetry and, thus, there is no symmetry rea- 
son to expect any fundamental relations among M's for different species unless additional 
assumptions are made. 

On the other hand, among relevant deformations to the dispersion relation the coef- 
ficient of the k'^ term should flow to the squared speed of light at the IR fixed point. (See 
Eq. (jlj).) This applies to all species, since relativity would not emerge in the IR, otherwise. 



Various aspects of Hofava-Lifshitz gravity were already explored in the literature 
example, cosmology based on Hofava-Lifshitz gravity was investigated in 



Some solutions in Hofava-Lifshitz gravit y w ere presented in 
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or 



extensions of the theory were proposed [15 
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IJ]. Possible 



The dispersion relation ([2]) leads to interesting cosmological consequences, such as gen- 
eration of scale-invariant cosmological perturbations [5|] and time-delays in Gamma-ray 



bursts [10|]. The purpose of this paper is to explore yet another cosmological consequence 
of the z = 3 dispersion relation. We shall point out that the radiation energy density in the 
UV epoch (T ^ M) is proportional to and, thus, decays faster than in the IR epoch 
(T <^ M) or in relativistic theories. This leads to intriguing cosmological consequences such 
as enhancement of baryon asymmetry and stochastic gravity waves. We shall also discuss 
current observational constrains on the dispersion relation. 

This paper is organized as follows. In Sec. [TT] cosmological impacts of dispersion relation 
([2]) are explained with some explicit examples. In Sec. IIIII observational constraints on the 
dispersion relation, in particular, the value of M is discussed. Sec. IIV]is devoted to summary. 



II. COSMOLOGICAL IMPACT OF DISPERSION RELATION 

We first derive the peculiar dependence of radiation energy density on the scale factor of 
the Universe. Suppose that the constituent particles satisfy the dispersion relation, 

where we set the speed of light c to be unity, for simplicity. If the typical momentum of the 
particle (~ temperature) is much lower than M, we recover the ordinary dispersion relation, 
uj'^ ~ fc^. On the other hand, if the momentum is much larger than M, the dispersion 
relation can be approximated to be 

k^ 

Assuming the relation ([5]), the energy density of the radiation, p^, decreases as 

Pr ~ <jjn oc a~^, (6) 

where a; is a typical energy of the particles, n the number density, and a represents the scale 
factor. This corresponds to the fluid with the equation of state = where is the 
pressure of the radiation. Here we have used a fact that the physical momentum k redshifts 
as A; oc due to the cosmic expansion, and the number density decreases as n oc . The 
peculiar dependence (EI) should be contrasted to the ordinary scaling, p^ oc a~^. Since the 
energy density of the radiation decreases more quickly than the ordinary case, the estimate 
on cosmological abundance of e.g. the baryon number, the stochastic gravity waves, etc. 
can be modified. Let us see this below. 
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A. Baryon asymmetry 



As an example we take up non-thermal baryon production from an inflaton decay. *^ 
Suppose that, after inflation ends, the inflaton cj) decays into radiation, which is heated up 
to a temperature higher than M. At the same time, the baryon asymmetry is assumed to 
3e generated from the inflaton decay. This is indeed the case of non-thermal leptogenesis 
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211 ]. As the Universe expands, the temperature will decrease. We assume that 



the ordinary dispersion, uj = k, is restored at a temperature below M. Then, the final 
baryon-to-entropy ratio can be estimated as 
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Reheating 



where 



denotes the inflaton energy density.*^ The final baryon asymmetry will be 

e T^Mp 

f^u _ 

T=M 



riB 
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(7) 



(8) 



M 

where Mp ~ 2.4 x 10^*^ GeV is the reduced Planck mass, denotes the decay rate of the 
inflaton, and e is the baryon asymmetry generated from the decay of one inflaton quanta. 
This result should be compared with the ordinary result. 



Up ^yv^Mp 

— ~ e . 

s 



Thus, we conclude that the baryon asymmetry is enhanced by 



A 



M 



(9) 



(10) 



compared to the case of the ordinary dispersion relation. Such an enhancement may make an 
inefficient baryogenesis scenario viable, which would result in too small baryon asymmetry, 
otherwise. 



B. Gravitational waves 



Following the similar arguments, it is also possible to enhance abundance of the relic 
gravity waves. Although the amplitude of the gravitational waves remain constant when 

We call a field as inflaton, which dominates the energy density of the Universe and decays into radiation. 
It does not affect our discussion whether or not the field actually induces the inflationary expansion. 
Note that the entropy density is proportional to for T < M, but this could be modified for T > M . 
The calculation here does not depend on the definition of the entropy for T > M. 



4 



the corresponding mode lies outside the horizon, once it enters the horizon its amphtude 
decreases inversely proportional to the scale factor. Since the timing at which a mode enters 
the horizon depends on the equation of state of the Universe, the thermal history of the 
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261 ] ■ The present 



Universe is imprinted in the gravitational wave spectrum 
spectrum of the stochastic gravitational wave background in terms of the density parameter 
rigviik) is given by 



^gw(^) — 



1 dp, 



gw 



2 A (prim) 
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where pco is the critical energy density at present, dpg^j/dlnk represents the energy density 
of the gravitational waves per logarithmic frequency, Hq is the present Hubble parameter, 
ain(A;) denotes the scale factor when the mode with comoving wave number k enters the 
horizon, and A^j^^^^\k) is the dimensionless power spectrum of the primordial gravitational 
waves, which is assumed to be scale invariant. From this equation we can see that the 
resulting spectrum scales as ilg^{k) oc A;(2-4p)/{i-p) jf ^j^g Universe expands as a cx when 
the mode k enters the horizon. Therefore we get 



^gw{k) oc < 



k-^ for p = 2/3 (w = 0), 
fco for p = 1/2 (w = 1/3), 
k^ for p = 1/3 {w = 1), 



(12) 



where w represents the equation of state of the Universe. The last line corresponds to the 
case where the radiation in the UV regime fills the Universe. Thus the mode with k > kM, 
where kj^ is the comoving Hubble scale at which the radiation transits from UV to IR 
regime, is enhanced by the factor A = k/kM, compared to the case of ordinary dispersion 
relation. The transition wavenumber kM is estimated as 



k 



M 



2.6 X 10'' Hz 



M 



(13) 



106.75 J \10^^ GeV^ 

where g^s denotes the relativistic effective degrees of freedom. As we will see in the next 
section, M is observationally bounded below as M > 10^^ GeV and hence the frequency 
( fT3|) lies outside the ran ge c overed by future space-based gravitational wave detectors, such 
as BBO and DECIGO 27|]. However, it may be possible to relax this lower bound as 
M > 10^ GeV as explained in the next section. In this case characteristic spectral shape 
may be observed by BBO/DECIGO. 
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FIG. 1: Present gravitational wave background spectrum. Solid (dashed) black line corresponds 
M = 10^^(10®) GeV. A mode enters the horizon in the matter dominated era, radiation dominated 
era with IR regime and radiation dominated era with UV regime as shown by MD, RD(IR) and 
RD(UV) on the top of the figure, respectively. Sensitivity of ultimate-DECIGO is also shown by 
blue dotted line. 

Fig. [1] shows the gravitational wave background spectrum. Solid (dashed) black line 
corresponds M = 10^^(10^) GeV. A mode enters the horizon in the matter dominated era, 
radiation dominated era with IR regime and radiation dominated era with UV regime as 
shown by MD, RD(IR) and RD(UV) on the top of the figure, respectively. Sensitivity 
of ultimate-DECIGO is also shown by blue dotted line. Here the tensor-to-scalar ratio is 
assumed to be 0.1.*^ It is seen that future gravitational wave detectors may have a chance 
to detect characteristic spectral shape of the gravitational waves if M happens to be close 
to the current observational bound. 

We have considered two simple examples above, but the enhancement can occur in other 
situations. For instance, if the dark matter (e.g. gravitino) is produced from the infiaton 



'^'^ Throughout this section we have assumed that the gravity waves satisfy the ordinary relativistic dispersion 
relation. This requires a scale Mgw, above which the dispersion relation for the gravity waves is modified 
like (HI, to be larger than the Hubble parameter during inflation. Since Mgw may be different from M, 
this is not a serious issue. If Mg^ is smaller than the Hubble parameter during inflation, the ampHtude of 
the gravity waves will be given by ~ Mg^/Mp instead of Hi^f/Mp, according to the discussion along with 
Ref. [5|]. Even if Mgw equals to M, the enhancement of the gravity waves still exists, although the detection 
of the stochastic gravity waves by DECIGO becomes slightly difficult due to smaller tensor-to-scalar ratio. 



6 



decay, its abundance will be similarly enhanced. Or, if preheating occurs soon after the 
inflaton begins to oscillate, significant amount of gravitational waves may be produced with 



peak frequency corresponding to the comoving Hubble scale after inflation 28|, |29|, 
These gravitational waves are also enhanced in a similar way. 
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III. OBSERVATIONAL CONSTRAINT ON THE DISPERSION RELATION 



In the previous section we have assumed that the radiation, including the standard-model 
particles, are subject to the dispersion relation ([S]) at a temperature higher than M. There 
are observational constraints on the possible deviation from the ordinary dispersion relation, 
using the arrival timing of the high-energy gamma rays coming from far distant sources 32| , 
such as gamma-ray bursts and active galactic nuclei. In fact, the dispersion relation (jl]) 
leads to the energy-dependent photon velocity. 



du 
dk 



ijj 



(14) 



where we have restored the speed of light c in the IR limit. In the IR limit {k <C M), we 
obtain 



3 k^ 



(15) 



If K is 0(1), the second term gives dominant contri 



speed. According to the MAGIC collaboration 

case reads*^ 



3ution to the variation of the propagation 



33| , the lower bound on the scale M in this 



M > lO^^GeV 



(16) 



for K = 0(1). The essence of the Hofava-Lifshitz quantum gravity is the appearance of the 
first term in the r.h.s. of (jl]), and so, the precise value of k may be different from order unity. 
In particular, if k, is negligibly small, it is fc^ term in Eq. f|T5|) that is subject to observational 
constraints. In this case the constraint on M becomes relaxed as*^ 



M > lO^GeV. 



(17) 



More or less similar bound was also found by the H.E.S.S. 3J| and Fermi 35] collaborations. 
Even tighter constraint may come from observations of ultra high energy cosmic rays with energy higher 
than 10^ GeV, which may push the bound as M > 10^^ GeV. 
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Note that, since the coefficient of fc'^-term in Eq. (fTSll is positive, it is not possible to attribute 
the observed delay in the arrival time of gamma-rays to the energy dependence in the speed 
of light .^'^ 

In order to see the enhancement of the relic gravity waves by the future experiments, 
M must be smaller than 10^ GeV or so, as we have seen in the previous section. This can 
be made compatible with the current observational constraint by assuming k ^ 1. On the 
other hand, the enhancement of the baryon asymmetry can be viable for a broader parameter 
space. 

IV. CONCLUSIONS 

We have shown that the non-trivial dispersion relation suggested by the Hofava-Lifshitz 
quantum gravity theory leads to a peculiar dependence of the radiation energy density on 
the scale factor of the Universe. Such a dependence leads to an enhancement of the baryon 
asymmetry, gravity waves, dark matter, and so on. We have demonstrated this enhancement 
in simple examples. 

These phenomenological aspects are not limited to the Horava-Lifshitz theory of gravity, 
but applied to any quantum gravity model in which the dispersion relation is modified at 
high energy. 
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